- Running a series of stratified analysis can be cumbersome
library(Rstrap)
## Load survey data
load("converted_set_details_2017-04-27.Rdata")
load("converted_length-frequency_data_2017-04-27.Rdata")
load("age-growth_data_2017-04-27.Rdata")
## 3NO Witch Flounder
witch <- strat.fun(setdet = setdet, data.series = "Campelen",
species = 890, survey.year = 1995:2012,
season = "spring", NAFOdiv = c("3N", "3O"))
## 3K Roughhead Grenadier (length disaggregated)
grenadier <- strat.fun(setdet = setdet, lf = lf, data.series = "Campelen",
program = "strat2 & strat1", species = 474, survey.year = 2012,
season = "fall", NAFOdiv = "3K", strat = 622:654, sex = "unsexed",
length.group = 0.5, group.by = "length")
## 2J3KL Greenland Halibut (length and age disaggregated)
turbot <- strat.fun(setdet = setdet, lf = lf, ag = ag, data.series = "Campelen",
program = "strat2 & strat1", species = 892, survey.year = 1998:2012,
season = "fall", NAFOdiv = c("2J", "3K", "3L"),
sex = c("male","female","unsexed"), length.group = 1,
group.by = "length & age")
## 3Ps Atlantic Cod (length and age disaggregated)
cod <- strat.fun(setdet = setdet, lf = lf, ag = ag, data.series = "Campelen",
program = "strat2 & strat1", species = 438, survey.year = 1998:2012,
season = "spring", NAFOdiv = "3P", sex = c("male","female","unsexed"),
length.group = 1, group.by = "length & age")
- This doesn’t include code for plotting these results